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Abstract. The partial structure factors (PSF) of molten Z11CI2 are known since 
1981 due to the neutron diffraction experiments by Biggin and Enderby. It is 
shown in this article that the set of PSF derived from this experiment is not 
consistent with new XRD data with much better statistical accuracy. A new 
set of PSF is derived which corrects this deficiency. The first diffraction peak is 
not dominated by the cation-cation-PSF and the set of PSF shows a remarkable 
similarity to the PSF of glassy GeC>2 and Si02 determined recently. The Levcsquc, 
Weis, Reatto potential inversion scheme has been used to interpret the PSF in 
terms of three dimensional structures. The thus derived angular correlations 
between near neighbor atoms show similarities with CuBr melt, reflecting the fact 
that both melts are based on tetrahedral structural units but also characteristic 

differences in the H h (ZnClZn) angle distribution, the angle joining adjacent 

tetrahedra. Further, a plausibility for the enhanced electrical conductivity of 
ZnCl2 at higher pressures is given, based on simulations with the equilibrated 
potential. 



PACS: 61.20.Qg, 61.20.Ja, 61.25 -f, 66.10.Ed 

1. Introduction 

ZnCi2 and other strongly structured divalent and trivalent chloride melts have been 
intensely studied in the past [1, 2, 3, 4, 5]. In ZnCl2 the tendency to form long ranged 
structures without periodicity is so pronounced, that it can be easily supercooled into 
a glass [6]. It has been found [7] |, that the ZnZn and the C1C1 next neighbor shell 
occur at the same distance. Such structures are especially difficult to describe by ionic 
models of the liquid structure. The reason for this difficulty is asymmetry between 
cation and anion charge. If such an ionic potential is applied, the structure tend to be 

I This reference is called I in the following 
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fmorite-like, i. e. the doubly charged ions further apart than the singly charged ions 
[8]. The inability of ionic models to reproduce the experimental structure has been 
taken as indication [6], that the bonding in ZnCl2 has a more covalent character. 

In the early eighties the partial pair distribution functions of ZnCl2 have been 
determined by Biggin and Enderby (I) with neutron diffraction and the isotopic 
substitution (NDIS) technique. Wilson and Madden, when forwarding their view that 
ZnCl2 can be described entirely by a polarizable ion model as opposed to a largely 
covalent description [9, 10] refer to this now almost twenty years old NDIS study to 
evaluate the success of their model to describe the ZnCl2 structure. A special accent 
has been given to the attribution of the first peak in momentum space to the ZnZn 
partial structure factor (PSF). However, the ZnZn-PSF, as it has a low weight in a 
total neutron structure factor, is the one which is most difficult to determine. It will 
be shown in this text, that the PSF determined by in (I) are not compatible with 
more recent neutron [11] and electromagnetic radiation scattering experiments [12] 
and an alternative set of PSF is determined which corrects for these deficiencies. The 
potential inversion scheme of Levesque, Weiss and Reatto [13] (LWR-scheme) is then 
applied to interpret the PSF in terms of three dimensional structures. 



2. Theoretical Summary 

The differential cross section of a liquid in a neutron or electromagnetic radiation 
scattering experiment can be expressed in terms of a total scattering function: 

s (n HQ) = Uj » " A +i (i) 
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(da_\( x ) I _ \r^N uc f2 

&'\Q) - HQ) + 1 = Uj % th + 1 (2) 

where (^) is the coherent differential cross section, Q — in/ A sin(6>), A the wavelength 
of the radiation and 9 the diffraction angle, hi the coherent scattering lengths [14], fi 
the X-ray form factors in the independent atom approximation [15], a e i the scattering 
cross section of the free electron, the stoichiometric coefficient of the atom i, and 
where the sums are extending over the number of distinct atoms N uc in the unit of 
composition, ZnCl2. These total structure functions are composed of partial structure 
factors PSF Sij\ 

SW*>=5>y(Q)*y(Q) (3) 

ij 

with 

™ ij (Q) = v -0^W 1 (4) 



where for the neutron case f(Q) has to be replaced by b and the weighting factors 
become independent of Q. The partial structure factors are related to the partial pair 
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distribution functions (PPDF) via Fourier sine transformation: 

r ■ ( 9ij - 1) = J Q ■ {an - 1) sin(Qr)dQ (5) 

with p uc the density per unit of composition. The PPDF gij describe the probability 
to find an atom of type j at distance r from an atom of type i relative to the bulk 
density. 

An alternative definition of the neutron total structure factor has been used by 
Biggin and Enderby: 

F(Q) = c 2 a b 2 a [S aa (Q) - l]+2c a c b b a b b [S ab (Q) - 1] +c 2 b b 2 b [S bb (Q) - 1](6) 

where c a , c b are the concentration of the atomic species, the S aa , S bb , S bb are partial 
structure factors differing from the Sij in Eq. 3 by a factor of v - It is noted, that in 
this definition the total neutron structure factor has the dimension of a cross section, 
while S^ n \Q) in Eq. 1 it is dimensionless. 

In order to generate three dimensional structures from the pair distribution 
functions, the potential inversion scheme of Levcsque, Weiss and Reatto [17] (LWR- 
scheme) has been applied. The idea of this method is based on the equation: 

-v(r) 



g(r) = cxp 



+ g(r) - 1 - c(r) + B(r, v) 



(7) 



kT 

relating the pair distribution function and the pair potential v(r), where c(r) the 
direct correlation function and B(r,v) the bridge function. Starting with a first guess 
of the bridge function, e. g. neglecting it completely, a first guess of the potential 
is calculated. A Monte Carlo (MC) or Molecular Dynamics (MD) simulation gives 
g(r)^ and c(r)( 1 ^ belonging to and thus B(r, v^ ). B(r,v^ ) is expected to be 
a better approximation for B(r, v) than completely neglecting it. Thus, substituting 
v) in equation 7 gives the LWR iteration formula 

v {n) j kT = v (n-l) + l n ( g («-l) /g^ x P)) 

+ cP l - 1] - cS exp) - g {n - 1] + g {exp) (8) 

The empirical potential Monte-Carlo (EPMC) scheme [16] is closely related to equation 
8, but considers the logarithmic term only. The complete form of Eq. 8 has been 
preferred here as this scheme shows faster convergence [17]. The applicability of the 
RLW scheme to polyatomic systems has been shown by [18]. An example of the 
application of this technique has been given recently in [19]. 



3. Deduction of the PPDF 

Biggin and Enderby (I) have determined the PSF and the PPDF of ZnCl 2 melt with 
neutron diffraction and isotopic substitution of the chlorine atom. This technique uses 
the fact that the neutron scattering power of the two chlorine isotopes for neutrons 
are significantly different. Thus, assuming the independence of the structure from the 
isotopes and the preparation involved, a set of three equations of the type given by 
Eq. 3 are obtained. 
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By coincidence the weighting factors Wij of the PSF in a neutron diffraction 
experiment on Zn 37 Cl2 and an electromagnetic radiation scattering experiment are 
the same, except for the small Q dependence of the ^Wij. Thus, S^ 37cl and S^nc\ 2 
are equal. Fig. 1 compares these quantities as well as S^ 37Cl and S^Ob wn i cn would 
result from the PSF given in (I). It is evident from this figure, that the PSF from (I) 
are incompatible with the X-ray result. 

Figure 1. Comparison the total neutron scattering function of Zn 37 Ci2 (I) and 

the electromagnetic radiation scattering function of ZnCb [12] 

F Zn 37 Ch (Q) = 0.3956 • 10- 28 m 2 [s£ i ) 37 - 1] (Eq. 6) (open diamonds, from 

I) is compared with 0.395610" 28 m 2 i(Q) (small crosses, from [12]) as well as 
F Zn 37 Ch (Q) (solid line) and 0.395610" 28 m 2 i(Q) (broken line), which would 
result from the PSF given in I. 

Wilson and Madden determined the PSF of the ZnCl2 system from a computer 
simulation using a polarizable ion model. In Fig. 2 the results of a similar comparison 
of the experimental electromagnetic radiation scattering function i(Q) with the result 
obtained for the simulated polarizable ion PSF via Eq. 3 and 4 are shown. Again, 
the incompatibility of the polarizable ion PSF with the electromagnetic radiation 
scattering experiment can be seen. The overestimation of the low-Q peaks is noted. 

Figure 2. Comparison of the experimental total X-ray structure factor i(Q) 
(solid line) with i(Q) calculated from the PSF resulting from the simulation of 
Wilson and Madden [9] (symbols) 

The structure of ZnCl2, as a two component system, is described by three PPDF/ 
PSF. A set of three independent equations 2 is hence necessary to determine the PSF. 
The first two are the hard X-ray i(Q)[12] and the neutron S(Q) of of ZnCk) in natural 
isotopic composition [11], the procedure to obtain the third independent information 
is described in the following. The ZnCl-PSF is dominated largely by the contributions 
of the first peak in real space at 2.3 A(cf. Fig. 4). This peak is separated from the 
rest of the structure in real space even in the total g^ x \r)[12\. The ZnCl-PSF of (I) 
was, hence, split by Fourier-filtering into the part originating from the first peak in 
real space and the remaining contribution: ("^SznCi — S^ n ' FP ^ + 5 , (™^ rem ) and S^ n ' FP ^ 
was removed from the ZnCl-PSF ,5'(™< rem ) shows oscillations exceeding the error level 
only for Q < 4A _1 and is set to zero for Q > 4 A -1 . s( FP ) has been determined 
from the hard X-ray total structure factor i(Q). A new sz n ci is then given by the 
contribution of the first peak from the hard X-rays and neutron derived s ( n > rem ); 
sznCi = s( x ' FP ) + s (". rem ). This gives the missing piece of information to determine 
all three PSF/PPDF. 

Relying on one partial structure factor from the NDIS experiment more than on 
the other may seem somewhat arbitrary, but two reasons justify the choice: First, the 
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SznCi is l ess affected by statisticaf errors then SznZn since its weighting factor is large 
and, second, s^nCi can ^ e ta k en horn the hard X-ray data. 

The resulting partial structure factors are shown in Fig. 3. It is noted, that the 
attribution of the first peak in Q space to sznZn, which seemed obvious from the old 
NDIS data is not reproduced here. Instead, the maximum at lA _1 is mainly a result 
of the sine wave in SznCl belonging to the first peak in r space. It is underlined, that 
the PSF in Fig. 3 are consistent with all available diffraction information and that a 
tree dimensional arrangement of atoms exist, which is consistent with these PSF (this 
is shown later in this paper). It is noted that the PSF of vitreous germanium-dioxide 
determined by Price et al [20] and the PSF of vitreous SiC>2 calculated from the ab- 
initio simulation by Sarntheim et al [21] show a remarkably similar behavior. It has 
been hypothesized [22], that the first diffraction peak in AX2 type glass formers is 
generally dominated by the cation-cation PSF. Together with the ZnCl2 melt, there 
are now three prominent glass formers of AX2 composition, where no indication for 
such a domination is found. 

Figure 3. Comparison of the simulated structure factor with the experiment at 
T=873K 

Left: Comparison of the experimental (cf. preceding section) partial structure 
factors (solid line) with the simulation (broken line). 

Right: Comparison of the experimental total structure functions i(Q) und F(Q) 
(data points), the recomposited experimental partial structure functions (solid 
line) and the simulated PSF (broken line). Solid and broken line are nearly 
indistinguishable. 



Figure 4. Partial pair distribution function of Z11CI2 from the experiment (Solid 
line) and the simulation (broken line). 

4. Simulation results 

In Fig. 3 and 4 the PSF and the PPDF deduced as described in the preceding section 
and the result of a MC simulation with the converged potentials from the LWR 
scheme is shown. Further, a comparison to the total neutron and X-ray structure 
functions is given (cf. Fig. 2). The total structure functions agrees with the simulation 
results almost perfectly except at the lowest Q-values. A slight overestimation of 
the periodicity in gznCl and gcici with respect to the experimental functions is to 
be remarked, transforming into differences of the height of the peak at 2A _1 in Q 
space. It is, however, possible that this is due to the imperfection of the experimental 
functions composed from quite different sources and the simulation correct for these 
imperfections to some extend, due to the restrictions imposed by the necessary space 
filling of a simulation box of given size. A snapshot of the simulated ZnCl2 melt is 
shown in Fig. 5. 
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Figure 5. Snapshot of simulated ZnCl2 melt. 

Dark spheres are Zn-Atoms and light spheres Cl-Atoms, adjacent atoms are 
connected with solid cylinders. Shown is neighborhood of the Zn-atom in the 
middle of figure up to the forth shell. 

In Fig. 6 the angular correlations between neighboring atoms resulting from the 
MC-simulation with the converged LWR potential are shown. A distance criterion 
has been applied to decide, whether atoms are neighboring, "bonded", or not. The 
cut-off radius for both ZnZn and C1C1 pairs, 4.8A, corresponds roughly to the first 
minima in the partial pair distribution functions. For the ZnCl PPDF the minimum 
between first and second maximum is very pronounced, and the choice of the cut-off 
radius, 2.8A, is evident. The — I — (ClZnCl) distribution, peaking at cos(#) = —1/3, 

correspond to the geometry of the ZnCLt tetrahedra, while the distribution 

shows a pattern characteristic for random dense packing (cf. Finney's model [23] for 
ZnCl2 glass). These two distribution functions are very similar to the corresponding 

— I — and of CuBr(l) determined by Pusztai and McGreevy [24] with Reverse 

Monte-Carlo (RMC) simulation. This reflects the fact that this liquid also shows a 
tetrahedral ordering - Cu(I) is isoelectronic to Zn(II). On the other hand, there is a 

marked difference in the H h distribution function, the angle connecting adjacent 

tetrahedra, which is much more pronounced in the ZnCl2 case. This makes plausible 
why a continuous random network is a meaningful description of the ZnCl2 [12] , but not 

the CuBr structure. Such a preference for certain H h angles restrict the possibilities 

to arrange the elemental tetrahedral units in a periodic lattice and makes it plausible, 
why the ZnCl2 but not CuBr melt can be easily supercooled into a glassy state. 

Figure 6. Angle distribution functions of adjacent atoms 

The pair potential being determined at one particular state point, it can be used 
in an NVT-Monte Carlo simulation to predict of the structure at state points which 
are not easily accessible to diffraction experiments, assuming the state independence of 
that effective pair potential. In Fig. 7 the changes in the pair distribution functions are 
shown, which are obtained from a Monte-Carlo simulation using the LWR potentials 
determined above, when increasing the density by 3.3% 6.6% or 10%. A ten percent 
density increase corresponds to an approximate pressure of 3000 bar (calculated from 
the compressibility at atmospheric pressure). The expected change in the X-ray total 
scattering function i(Q) is shown in Fig. 7. These are obviously small, the most 
important change is the decrease in intensity at low-Q compatible with a decrease in 
compressibility to be expected at higher pressures. The ZnCl coordination number 
has a slightly increasing tendency raising from 3.6 to 3.7 at 10% compression. 

There has been recent interest to understand the behavior of the electric 
conductivity in ZnCi2 melt and strongly structured molten salts in general [5, 25, 26]. 
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Figure 7. Isothermal compression of Z11CI2 melt at 873 K. 

Left: Electromagnetic radiation total structure factor i(Q) at atmospheric 
pressure (solid line), 3.3%, 6.6% and 10% increased density (dotted, dashed dotted, 
dashed line) 

Right top: ZnCl-PPDF at atmospheric pressure (triangles), 3. 3% (crosses), 6.6% 
(stars) and 10% (diamonds) increased density. 

Right bottom: Change in the ZnCl-PPDF upon 10% compression 

Contrary to the usual trend, in ZnCl2 melt the conductivity increases with increasing 
pressure [26]. A Monte-Carlo simulation does not allow, of course, the calculation of 
dynamic correlation functions like the velocity correlation function directly related 
to the conductivity. However, it is assumed [5] that a velocity determining step 
in the ionic conduction in ZnCl2 is the detachment of the chlorine ion from the 
network. While the changes in the pair distribution in functions are small, the depth 
of the minimum between the first and second ZnCl coordination shell decreases. It 
is probable, that small number of atoms at intermediate distance between the first 
and second neighbor shell play a crucial role for the conductivity and correspond to 
detached atom in interstitial sites (A similar argument can be found in [27, 28]). Thus 
a plausible explanation of larger conductivity is that these interstitial sites become 
accessible for more ions and thus the number of ions accessible to charge transport 
increases. 

5. Conclusion 

Using the results of modern diffraction experiments improved partial structure factors 
and pair distribution functions for ZnCl2 melt are obtained. The first peak in Q-space 
at lA _1 has been traditionally attributed to the ZnZn partial structure factor. This 
assumption is not in agreement with the hard X-ray diffraction experiment. A set of 
PPDF consistent with all available diffraction experiments is presented here, which 
suggests that this peak is mainly due to the ZnCl partial. It is also proved in this 
work, that this set of PSF can be mapped to a physical arrangement of atoms of the 
experimental density. It is shown, that the simulation result of Wilson and Madden, 
although compatible with the old NDIS data of Biggin et al within its large statistical 
errors, shows significant deviations from modern hard X-ray data. 

In order to interpret the partial pair distribution function it terms of a three 
dimensional structure the Levesque, Weis, Reatto scheme has been used to deduce 
an empirical potential for this system. MC simulations with this potential are used 
to derive the angular correlations between neighboring atoms. It is interesting to 
compare the angle distribution functions of the ZnCl2 melt, which can be easily 
supercooled into a glassy state to those of the melt of the fast ion conductor CuBr. The 
structure of both systems is characterized by tetrahedral structural units, translating 

in a similar — I — distribution. On the other hand, the H \- distribution , the angle 

joining adjacent tetrahedra, is peaked in ZnCl2 and broad in CuBr, explaining why a 
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continuous random network is useful characterization of only the ZnCl 2 melt. 

Besides of the capability to produce three dimensional models of the ZnCh in 
excellent agreement with the available diffraction information, the potentials obtained 
by the LWR scheme can be used to predict the structure at state points difficult 
to access by diffraction experiments. Thus, the structure of the melt at elevated 
pressures has been obtained. It is found, that under these conditions, the ZnCl 
coordination number increases slightly, while the overall structural changes are small. 
The increasing ionic conductivity of ZnCl2 melt at higher pressures can be explained 
by the increased number of ions capable to occupy interstitial sites in the network. 
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